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Abstract 

We examine solitary waves in classical ferromagnetic Heisenberg chains with an uni- 
axial anisotropy and a parallel magnetic field in a continuum approach. The boundary 
conditions commonly used are generalized to nonlinear spin wave states, which themselves 
turn out to be stable only for an anisotropy of the easy-plane type. In this case we obtain 
two different branches of one-soliton-solutions which can be mapped onto each other by 
a formal time inversion. Moreover, they show some remarkable similarity to dark solitons 
of the Nonlinear Schrodinger equation. Numerical simulations for the discrete Heisenberg 
chain show that these solitary waves are highly, but not absolutely stable under interaction 
with linear excitations and as well under scattering with each other. The possible signifi- 
cance of these solitary waves in a phenomenological theory of one-dimensional magnets is 
briefly addressed. 

PACS numbers: 75.10.Hk, 03.20.-fi, 03.40.Kf 

1 The model 

The discrete classical ferromagnetic Heisenberg chain with a local uniaxial anisotropy and a 
magnetic field parallel to this axis is given by 



SnSn+l + hS^ + — (5^)^ 



(1) 



The classical spins are unit vectors Sn = (sini^^ cos (/?„, sinT?^ sint/?^, cos h the magnetic 
field and a the anisotropy parameter. If we choose the infinite chain to lie in the ^-direction 
and an appropriate length unit with lattice spacing a = 1, the continuum approximation of 
dl]) in lowest order reads 



Here the classical spins are described by the canonical conjugate fields p{x,t) = cos(i?(x, t)), 
q{x,t) = ip{x,t). The equations of motion are 



]- o2 P_ 

dtp = (1 - p^)dlq - 2p{d,p)d,q . (4) 



dtq = 2^1p - -pr^-^{dxpf - p{dxqf - h - ap , (3) 
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For a < and | — | < 1 the solution of lowest energy is given by constant q and p = —— (easy- 
plane or, for finite h, easy-cone model), for q > by |p| = 1 (easy-axis model). These field 
configurations can be used as boundary conditions for the above equations. In the literature 
only these two cases appear to be considered, for a review see ||l|, |2|. 

Takhtajan Q and Fogedby |^] have shown that the isotropic model (a = 0) allows for a 
Lax representation and is solvable by the inverse scattering method. Moreover, this model is 
integrable in the sense that it has an infinite series of independent conserved quantities. These 
results have been extended by Sklyanin to the case of a general biaxial anisotropy. The 
relationship of the model (|2|) to the Nonlinear Schrodinger equation has also been established 
for the isotropic |^, 0] and the anisotropic case |^, ^J. However, all these results have been 
obtained using boundary conditions of the above type. In the present work we examine the 
model for more general boundary conditions to be specified below. 

If the spin configuratiuon is the same at both boundaries, the system has (at least) two 
well-known conserved quantities apart from the energy, namely the total momentum 

P = J dxpdxq (5) 

and the z-component of the total angular momentum or magnetic moment 

M = y" dxp. (6) 

P is the generator of translations and M is the generator of uniform rotations of the spins. 
The Poisson brackets of TC, P, M with each other vanish since surface terms do not contribute 
under the above condition. 



2 Nonlinear spin waves and linear magnons 

From the above equations one easily finds the solutions 

p = u , q = Qo + kx — Lot (7) 

with constants u, k, qQ, \u\ < 1, and uj given by 

uj = (^k^ + a^u + h. (8) 

Here the spin field describes a propagating wave characterized by the parameters u and k. 
(|7|) is an exact solution of the nonlinear equations of motion and will therefore be called a 
nonlinear spin wave. It has a homogenous energy density obtained from (|^ as 

e = ^i^l-u^^k^-hu- . (9) 

Note that the spin configurations mentioned in the previous section are included as special 
cases. The time-dependent field configuration can be used as a generalization of the 
boundary conditions given above. In the next section we will obtain solitary waves obeying 
such boundary conditions, i. e. having the asymptotic form of a spin wave for |x| — > oo. As 
the spin waves in general do not give absolute minima of the Hamiltonian their stability is 
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to be questioned. Therefore let us briefly discuss the stabihty of the homogenous spin wave 
solution ^ if the system is forced to have the asymptotic structure 



lim p = u , lim dxQ = k . (10) 

>oo >oo 

The ansatz 

p{x,t) = u + r]{x,t) , q{x,t) = qo + kx — Ldt + ^{x,t) (11) 

with T], ^ vanishing for |x| ^ oo describes small fluctuations around the spin wave solution 
respecting the boundary condition and the conservation of P and M. Introducing the Fourier 
transforms 

r,{x,t) = dlf,{l,t)e''\ (12) 

e(x,i) = -^Jdliil,t)e^'- (13) 

this leads to the restriction fj^O, t) = 0. Inserting these expressions in the equations of motion 
and linearizing in r], ^ gives 

dtfi = -2iuklfi - f {l - u'^^ I , (14) 

dti = I ^ ^ - - a\ fj - 2iukll . (15) 
\ 1 — / 

The general solution of these equations has the form fj{l,t) = r7o(/)e~*'^'*, ^(/,t) = ^0(0^"*"^'* 
with 

uJi = 2ukl ± ^l* + l'^{l-u'^)[-a-k'^) . (16) 

The nonlinear spin wave solution is stable if uJi is real for every Fourier mode I. This is the 
case either for \u\ = 1, or for 

- a - /c^ > (17) 



with general u. The relation (|1^ can only be fullfilled for flnite A: if a < 0, and stable 
exact spin wave solutions ^ are restricted to such type of Hamiltonian. If ( [TTD is valid 
eq. ([l^ ) is the spectrum of linear excitations above the nonlinear spin wave. We shall call 
these excitations magnons. The extremal group velocities within the two magnon dispersion 
branches are achieved at / = with 

^\ = 2uk ± J{l-u'^){-a-k'^) . (18) 
dl J ' - ^ 



1=0 



The upper sign corresponds to a minimum, the lower to a maximum of dcoi/dl on the particular 
branch. 

3 Solitary waves 

As the nonlinear spin waves are stable uniform solutions under the condition (0) we use them 
as a reference state to construct localized solutions. We are looking for solitary solutions which 
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have the form of a spin wave for |x| — > oo determined by parameters u, k. 
Inserting the ansatz 

p{x,t) = p{x — vt) , q{x , t) = ^It + q{x — vt) (19) 
and performing two formal integrations one ends up with (cf. Fogedby Q) 



dx ) 



dq r — p 

dx 1 — p 

2 



^ , (20) 



(l + - 2rp) + 2 {h + n) p [p"^ - 1^ - s [p^ - l) + "P^ {p^ - l) > (21) 



where r, s are integration constants, v the velocity of the solitary wave and an internal 
frequency. Denoting the r.h.s. of (|2l|) by F{p) we conclude that this quantity should be 
nonnegative for certain values of p between (—1) and 1. A real single root of F{p) corresponds 
to an extremal value of p{x), while a real double root p = u leads to a fixed point since one 
can easily show that in this case also the second and successively all higher derivatives of p 
vanish for p = u. For our solution to be localized and asymptotically in the state p = u we 
need the latter case and make the general ansatz 

(^) ' = " + " ""^^ + " ""^^ ■ ^^^^ 
To fullfill the additional boundary condition for q we must also have 



r — u 

|a;|^oo i — U 



lim dxq = v- ^ = k . (23) 



A comparision of eqs. (poD, ( pij ) and (p2D, ( p3| ) leads to 

vr = k(l - u^'j + uv , (24) 

.2 



s = + 2kuv + k^ - 'iu^k^ - au'' (25) 
and 

fi = -2u (-a - fe^) +2k{v- 2uk) , (26) 

7 = - (v - 2uA;)2 + (-a - A;^) (l - n^) , (27) 

where the parameters v and are related by 

vk= (k^ + a)u + h + n. (28) 



The additional roots u± of F{p) are given by 



/3 /32 7 

U+ — U = ± \ / - — ?r 

2a V 4a2 « 



ti/c^ vk , -a-{v-ukY , , 

-M + ±W ^ '-{-a-k^) (29 

a a V a 
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3.1 The case a <0 



The possible behavior of F{p) for a < and real u± is shown in figure |T|. For a localized 
solution we need n_(_ > n > u_ (a), while otherwise, e. g. case (b), one obtains a periodic 
wave train oscillating between n_|_ and U-. If a < a necessary and sufficient condition for 
'W+ > w > ti_ is 7 > 0, which also ensures the reality of u±. Thus we have 

7 = - (t, - 2ukf + (-a - (l - n^) > (30) 

as a restriction to the soliton parameter v. Moreover, the desired solution is only meaningful 
if it holds \p{x)\ < 1 for all x, i. e. < 1. Using a < one can derive from the above 

equations that < 1 is equivalent to {v — k{l + u))'^ > and U- > —1 is equivalent to 
(v - k{l - ■u))2 > 0. Thus, for every choice of the soliton parameter v compatible with (pQ) 
we have two solitary solutions which can be obtained explicitly by elementary integration of 
eqs. (HI), (H: 

P±i^,t) = n + - ^ J±^±-^ ^ (31) 

(7 + ^ {u± - u)j cosh {^yJ{x -vt- Xq)) — § {u± - u) 

q±{x, t) = qo + kx — wt 

_i f V + k (1 — u) u± — u , /a/7/ n\\ 
— tan — tanh - — (x — vt — xq) ] 

V V7 i + ^± V 2 ^ "Vy 

_i f V — k (I + u) u± — u /a/7/ 
-tan- ^ — -^—-tanh^—ix-vt-xo))j 

(32) 

The double sign corresponds to a different choice of the integration constant p{xQ,0) = u±, 
and go, xq are further constants. The two solutions are related via a formal time inversion, 
i. e. the mapping {x,t) 1— > (x,— t), {q,p) 1-^ iq,—p)i {h,a) ^ {—h,a), {k,u) ^ {k,—u), 
{v,Q) ^ (— u,— r^). As a consequence we have (/3, 7) 1— > (— /?,7), u± ^ —u^, and the (+)- 
and the (— )- solution interchange. Note also that the Hamiltonian and the equations of 
motion keep their form under such operation. 

We have found the general solitary wave of the form (|l9|) obeying the boundary conditions 
([lO|) . The solutions for p are pulse solitons with amplitude {u± — u) and width 1 / -^7 which 
are parametrized by the velocity v; a typical example is given in figure |2|. Prom (|3l|), ( |32[ ) one 
explitly sees the structure of a spin wave (0) for |x| » l/\/7- In this sense these solutions may 
be called dark solitons. The phase velocity of the spin wave and the soliton velocity differ by 
n/k (cf. (P), (^)), which can alternatively be used as soliton parameter. Note also that the 
range of admissible velocities given by (|3^) is centered around the formal group velocity of 
the spin wave. Moreover, as seen from eq. (^) the maximum (minimum) soliton velocity is 
given by the minimum (maximum) group velocity of the linear magnons on the upper (lower) 
magnon branch. 

The above results are in agreement with earlier work by other authors, beginning with Akhiezer 
and Borovik |10|, who took the absolute ground state of the model as boundary condition. 



For \u\ < 1 the inequality (|30| ) can only be fullfilled if 

- a - /c^ > . (33) 
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This is identical with the condition ( p^ ) for the stabihty of the spin waves. Thus, the above 
pair of one-sohton~solutions exists exactly in the region of boundary parameters u, k where 
the asymptotic structure of the solutions is found to be stable. 

If V takes its extremal values, i. e. 7 = 0, and assuming /? > 0, the (— )-soliton becomes a 
simple spin wave, while the (+)-soliton gets an algebraic structure: 



p+ix,t) 
q+{x,t) 



1 + -^{x -vt- Xof 
qo + kx — ujt 

' (3 {v + k{l — u)) {x — vt — xq) 



(34) 



tan 



tan 



—a + (3 — au 2 
I3{v-k{l + u)) {x-vt 



P + au 2 J ^^^^ 

If /3 < for an extremal value of v the above results hold vice versa. An algebraic form of 



-a 



a soliton in such limiting cases has also been obtained by Ivanov et al. [11 1 considering the 
special case u = k = 0. 

The energy density of the solitary waves is given by 

e± = —a {p± — uf + {—k {v — uk) — h — au) {p± — u) , (36) 

where we have subtracted the homogenous energy density of the underlying spin wave. From 
the explicit solutions (|3T[), (32) we calculate the total energy and the quantities P, M: 

2h f , f (3 \ . TT 



E± 



2^/7- 



-a 



tan 



(37) 



dx {p± {dxq±) - uk) 
(1+u)- 



tan 



-1 



2 



7 



TT 



tan 



,/^{v + k{l-u))^ 
n) + 7 



1(1 



^{v-k{l + u)) 



lb —sign {v + k{l 



± -Sign [V 




dx {p± — u) 



tan 



IT 

± - 

2 



(38) 



(39) 



-a \ \2^/^aJ, 

pi), (|39|) we have also subtracted the homogenous contributions of the nonlinear 



In eqs. 

spin wave from the integrands in (^), (^). For vanishing magnetic field and a given velocity 
both solutions have the same energy. In particular, for h = 0, (3 > the limiting (+)- 
soliton (|3^ , (|35|) has the same energy as the (— )-soliton, i. e. the homogenous spin wave 



Eq. (37), (38) provide a parametric representation of the dispersion law E{P) in terms of the 
soliton velocity. In figure ^ we give an example for the function E{P) for both soliton types. 



3.2 The case a > 

For a > and 1^1 < 1 we have shown that the spin wave solutions are unstable. Moreover, as 
the integration of eq. (p2|) is formally the same as for a < 0, we see from (|3l|) that we need 
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7 > for a localized solution. Because of ( p3| ) this cannot be achieved for nonnegative a. 
Thus no localized solution of the form (|l9|) exists for |tt| < 1 

The case \u\ = 1 is particular. The corresponding spin wave realizes the absolute minimum of 
H and the wavenumber k becomes irrelevant. For such boundary condition well-known solitons 
exist, which can be derived similarly by elementary integration (Kosevich et al. ||l^, Long 
and Bishop These solitons are characterized by two independent parameters v and $7 

because the relation (28) is not valid if the boundary condition ( p^ ) is not imposed. Therefore 



the soliton can alternatively be described by two independent conserved quantities, the total 
momentum and the total angular momentum, and the exact dipersion law E = E{P, M) can 



be obtained fl^, 14 1. These two-parametric solutions are also valid for a < and suffciently 



large {h-\-Vt), and under further conditions algebraic solitons arise as limiting cases |]TT[. 

4 Stability and scattering of solitons 

We now want to establish the stabitlity of the solitary waves under interaction with magnons 
and under scattering with each other. The stability is clear if one can prove the integrability 
of the model for the generalized boundary conditions used here. The solitary waves have some 
remarkable similarity to dark solitons occuring as solutions to the Nonlinear Schrodinger 
equation with repulsive interaction |15|. The latter evolution equation is well-known to be 



integrable ||T^. Moreover, Nakamura and Sasada have proposed a gauge transformation 
mapping the Nonlinear Schrodinger equation with attractive (repulsive) interaction on the 
Heisenberg chain with easy-axis (easy-plane) anisotropy, where the boundary conditions to 
the spin models correspond to the absolute ground state. Unfortunately, for the easy-plane 
case which is of interest here, this work suffers from some errors pointed out by other authors 



1 17, Thus, we shall use other means to examine the question of stability. 
In order to provide a brief demonstration of the physical significance of the solitary waves 
presented before and also to study discreteness effects to these continuum solutions we have 
carried out numerical simulation for the discrete model (|l]). The spin dynamics is given by 
the following Landau-Lifshitz-equation: 

dtSn = 5„ X [Sn-i + Sn+i) + /i (5„ X e,) +a(Sn- e,) (Sn X e,) (40) 

with S"^ = 1 and being the unit vector in z-direction. This equation may also be written in 
an explicitly Hamiltonian form, e. g. using the canonical conjugate variables pn = cosi?„, qn = 
ipn- It is seen easily that the discretized version of the nonlinear spin wave, i. e. pn = u, qn = 
kn + ujt + const., is an exact solution of the dicrete model ]l9[. The natural time unit of the 
system is chosen to be unity by the definition (|l]) . We numerically integrated the equations of 
motion (|40| ) on lattices of up to 5000 sites using a Runge-Kutta scheme of fourth order with 
a stepwidth of 0.01 . . . 0.05. The boundary conditions were implemented by rotating the spins 
on the terminal lattice sites with a constant frequency uj given by (|8|) while the z-component 
has the constant value u. 

If one uses the one-soliton solution of the continuum model evaluated at the discrete lattice 
sites as initial data for a simulation, magnons are emitted from the soliton during the first 
few hundred time units. This is an discreteness effect since the continuum solitary wave can 
only be an approximate solution to the discrete model. After this process the soliton moves 
with constant velocity over the lattice, with no measurable sign of instability over several 
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thousand time units. The separation of magnons and sohtons corresponds to the result in 
the continuum model that the admissible soliton velocities and magnon group velocities are 
strictly separated. In figure ^ we show an example where this effect is comparatively strong 
for a certain choice of parameters. Moreover, small magnons are also radiated from the edges 
of the system during the simulation. 

Since the spin configuration outside a soliton decreases rapidly to the underlying nonlinear spin 
wave, different one-soliton solutions can be matched and scattering experiments can be carried 
out. If one uses appropriately matched continuum solitons as initial data, the two solitons 
can apparently pass each other keeping perfectly their shape and identity. More precisely, we 
found no effect of instability on the background of magnons that stem from the adjustment 
to the discrete system in the beginning of the simulation or are radiated from the edges. This 
holds for all values of the system parameters and energy of both (+)- and (— )-solitons. A 
typical example is shown in figure ^. 

For more precise measurements the magnon background has to be reduced. This can be done 
by using the soliton shape after the adjustment to the lattice, i. e. the emission of magnons, as 
initial data for further simulations. Different spin configurations obtained by this procedure 
can be matched for a scattering experiment. This method leads to much purer initial conditions 
and works particularly well in the case h = 0, k = 0, where the spin configuration outside the 
solitons is spatially constant. In order to prevent the small magnons from the edges disturbing 
the scattering process one can choose a sufficiently large system. In figure ^ we show a typical 
scattering process. On the natural scale of the problem the two solitons apparently pass each 
other without any sign of instability. If the scale of the vertical axis is magnified by a factor 
of about 10^ one clearly sees that the solitons do not perfectly keep their shape but magnons 
with amplitudes of order 10~^ are emitted as a result of the scattering process. We have 
repeated this experiment with different stepwidths of the Runge-Kutta scheme and confirmed 
that this effect is not a numerical artifact. 

This observation strongly indicates the nonintegrability of the discrete model in the strict 

mathematical sense. Of course, this does not imply any statement about the integrability of 

the continuum model under the generalized boundary conditions. In particular not, since the 

discrete model shows properties very near to those of integrable models. 

In summary, we have demonstrated by numerical simulations that our solitary waves are 

physically significant objects showing a highly distinct but not absolute stability in their time 

evolution. 

5 Conclusions 

In this work we have examined the classical Heisenberg chain with uniaxial anisotropy under 
generalized boundary conditions. For an anisotropy of easy-plane type two branches of soli- 
tons are found on the background of a nonlinear spin wave that is used as boundary condition. 
The parameter area of stability of the nonlinear spin wave solution coincides with the exis- 
tence of the solitary waves. For the easy-axis and the isotropic model we have shown that no 
further one-soliton solutions exist apart from those already known. Moreover, in these cases 
the spin wave solutions are found to be unstable. This fact has apparently not been realized 
yet, since the spin wave solutions (^) have been discussed for such systems by various authors 
without mentioning their instability [||, ^, ^ . 
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For the case a < we have examined the stabihty of the continuum sohtary waves by numeri- 
cal simulations of the discrete Heisenberg chain. The solitary waves turn out to be highly but 
not absolutely stable in their time evolution and especially under scattering with each other. 
Thus the discrete model is shown to have properties very near to those typical for integrable 
models, but small effects indicating nonintegrability are clearly observed. From these numer- 
ical results and also from the similarity of the solitary waves obtained here to dark solitons 
of the Nonlinear Schrodinger equation one can strongly conjecture the integrability of the 
continuum model. 

The spin wave state could serve as a phenomenological model for a one-dimensional easy- 
plane magnet at finite (not necessarily very low) temperature. In thermal equilibrium the 
spatially averaged energy density and magnetization have certain values determining the pa- 
rameters u, k. However, the long-range ordering in this one-dimensional system must be 
destroyed at finite temperature by thermal fluctuations |2l|. The linear and localized nonlin- 
ear excitations obtained in this work provide a possible mechanism for this effect. 
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Figure 1: The typical behaviour of the polynomial F{p). 
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Figure 2: Spin field configuration of a (-|-)— soliton with system parameters a = —0.5, h = 
0.1, u = 0.4, k = 0.2 and soliton parameters v = 0.4, xq = 0, qo = as a function of x at time 
t = 0. 
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Figure 3: Dispersion law E{P) for particular values of the system parameters. The branches 
of (+)- and (— )-solutions meet in the origin. The dashed line is a guide to the eye. 
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Figure 4: Time evolution with a continuum solitary wave as initial data (small graph, (— )- 
solution) on a lattice of 1000 sites. For this choice of parameters comparatively large magnons 
are emitted in the first few hundred time units; only the component is shown. 
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Figure 5: Soliton scattering in a system with a = —0.5, h = 0.2, u = 0.5, k = 0.25; only the 
component is shown. 
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Figure 6: A scattering experiment with very pure initial conditions. In the insets the z- 
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